Entanglement entropy in the Lipkin-Meshkov-Glick model 
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We analyze the entanglement entropy in the Lipkin-Meshkov-Glick model, which describes mu- 
tually interacting spins half embedded in a magnetic field. This entropy displays a singularity at 
the critical point that we study as a function of the interaction anisotropy, the magnetic field, and 
the system size. Results emerging from our analysis are surprisingly similar to those found for the 
one-dimensional XY chain. 
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Within the last few years, entanglement properties of 
spin systems have attracted much attention. As initially 
shown in onc-dimensional (ID) spin chains 0, 0, 0, Of, 
observables measuring this genuine quantum mechanical 
feature are strongly affected by the existence of a quan- 
tum phase transition. For instance, the so-called concur- 
rence 5] that quantifies the two-spin entanglement dis- 
plays some nontrivial universal scaling properties. Sim- 
ilarly, the von Neumann entropy which rather charac- 
terizes the bipartite entanglement between any two sub- 
systems scales logarithmically with the typical size L of 
these subsystems at the quantum critical point, with a 
prefactor given by the central charge of the correspond- 
ing theory 0, 0, la • Note that the role played by the 
boundaries in these conformal invariant systems has been 
only recently elucidated [a] . Apart from ID systems, very 
few models have been studied so far [3, LLI1 [llj, LU| , 
either due to the absence of exact solution or to a dif- 
ficult numerical treatment. In this context, the Lipkin- 
Meshkov-Glick (LMG) model O EE El 

discussed here 

has drawn much attention since it allows for very effi- 
cient numerical treatment as well as analytical calcula- 
tions. Introduced by Lipkin, Meshkov and Glick in Nu- 
clear Physics, this model has been the subject of inten- 
sive studies during the last two decades because of its 
relevance for quantum tunneling of bosons between two 
levels. It is thus of prime interest to describe in particu- 
lar the Josephson effect in two-mode Bose-Einstein con- 
densates 0, 0. The entanglement properties of this 
model have been already discussed through the concur- 
rence, which exhibits a cusp-like behavior at the criti- 
cal point 0, Eol pll Eij l as well as interesting dynam- 
ical properties |2j|]. Note that similar results have also 
been obtained in the Dicke model [2J, |25| which can be 
mapped onto the LMG model in some cases [2(|, or in 
the reduced BCS model H3- 

In this letter, we analyze the von Neumann entropy 
computed from the ground state of the LMG model. We 
show that, at the critical point, it behaves logarithmi- 
cally with the size of the blocks L used in the bipartite 
decomposition of the density matrix with a prefactor that 



depends on the anisotropy parameter tuning the univer- 
sality class. We also discuss the dependence of the en- 
tropy with the magnetic field and stress the analogy with 
ID systems. 

The LMG model is defined by the Hamiltonian 



A 



(i) 



i<j 



where a\ is the Pauli matrix at position k in the direction 
a, and N the total number of spins. This Hamiltonian 
describes a set of spins half located at the vertices of a 
iV-dimensional simplex (complete graph) interacting via 
a ferromagnetic coupling (A > 0) in the xy-spin plane, 7 
being an anisotropy parameter and h an external mag- 
netic field applied along the z direction. The Hamiltonian 
l|T)l can be written in terms of the total spin operators 



H 



A(i + 7) (s 2 -^ 



N/2) - 2hS z 



~(l-7)(fi 



S 2 



(2) 



In the following, for simplicity we set A = 1, and since 
the spectrum of H is even under the transformation h «-> 
—h I23, we restrict our analysis to the region h > 0. 
Furthermore, we only consider the maximum spin sector 
S = N/2, to which the ground state is known to belong. 
A convenient basis of this subspace is spanned by the so- 
called Dicke states \N/2,M), which are fully symmetric 
under the permutation group and are eigenstates of S 2 
and S z , with eigenvalues iV(iV+2)/4 and M, respectively. 

In order to study the entanglement of the ground state, 
we need to define an appropriat e fig ure of merit. Follow- 
ing the ideas presented in [3, |j, |lfj, , we consider the 
von Neumann entropy associated to the ground-state re- 
duced density matrix pl.n of a block of size L out of the 
total N spins, Sl,n = — Tr (pl.n ^og 2 Pl.n) and analyze 
its behavior as L is changed, both keeping finite or 
sending it to infinity. Note that since the ground state 
reduced density matrix is spanned by the set of (L + 1) 
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entanglement entropy 




FIG. 1: Entanglement entropy for N 
function of h and 7. 
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Dicke states, the entropy of entanglement obeys the con- 
straint Sl,n < log 2 (i + 1) for all L and N, where the 
upper bound corresponds to the entropy of the maxi- 
mally mixed state pl,n = ^/(L + 1) in the Dicke basis. 
This argument implies that entanglement, as measured 
by the Von Neumann entropy, cannot grow faster than 
the typical logarithmic scaling observed in ID quantum 
spin chains [JQ. 

In order to study the different entanglement regimes, 
we compute the entropy in the plane spanned by 7 and 
h. This numerical computation can be laid out analyti- 
cally taking advantage of the Hamiltonian symmetries to 
reduce the complexity of the numerical task to a poly- 
nomial growth in TV. Results are displayed in Fig. \I\ior 
N = 500 and L = 125. For 7 / 1, one clearly observes 
an anomaly at the critical point h = 1, whereas the en- 
tropy goes to zero at large h since the ground state is 
then a fully polarized state in the field direction. In the 
zero field limit, the entropy saturates when the size of the 
system increases and goes to S^n = 1 f° r 7 = 0, where 
the ground state approaches a GHZ-like (cat-like) state 
as in the Ising spin chain |3|,|J]. By contrast, for 7=1, 
the entropy increases with the size of the system in the 
region < h < 1 and jumps to zero at h — 1 as we shall 
now discuss. 

In the isotropic case (7 = 1), it is possible to analyti- 
cally compute the entropy of entanglement since, at this 
point, H is diagonal in the Dicke basis. The ground-state 
energy is given by E (h, 7=1) = + j^M 2 - 2hM, 
with 



M 



I(hN/2), for < h< 1 
N/2, for h > 1 



(3) 



and the corresponding eigenvector is simply \N/2,M). 
Here, I(x) denotes the round value of x. 



To calculate the entropy, it is convenient to introduce 
the number n of spins "up" so that M — n — N/2, and to 
write this state in a bipartite form. Indeed, since Dicke 
states are completely symmetric under any permutation 
of sites, it is straightforward to see that the ground state 
can be written as a sum of byproducts of Dicke states: 



\N/2,n-N/2) 



L 

^p\ ,2 \L/2,l- L/2) ® 
\(N-L)/2,n-l-(N- L)/2) 



(4) 



where the partition is made between two blocks of size L 
and (N — L) and 



Pi 



N - L 
n — I 



(5) 



The ground-state entropy is then simply given by 

SL,N{h,^) = -Ya=oPi\°&2Pi- In tnc umit N , L > 
1, the hypergeometric distribution of the pi can be 



recast 
1 



into a Gaussian distribution 



Pi 



pi 



a _ 



2lTO 



exp 



of mean value I 



n-fi and variance 



n(N — n) 



(N — L)L 

iV 5 ' 



where we have retained the subleading term in (AT 
explicitly preserve the symmetry Sl,n = Sn 
entropy then reads 



L,N 



(6) 

L) to 
The 



f°° 1 
- / dl pf log 2 pf = - (log 2 e + log 2 2tt + log 2 a 2 ) , 
J -00 z 

(7) 

and only depends on its variance as expected for a Gaus- 
sian distribution. Let us mention that this result has 
been recently obtained in the context of the ferromag- 
netic Heisenberg chain j^. Of course, for h > 1, the en- 
tanglement entropy is exactly zero since the ground state 
is, in this case, fully polarized in the magnetic field direc- 
tion (n = N). For h <E [0, 1[ and in the limit N,L > 1, 
Eqs. ©, © and Q) lead to 



1 



S L .N(h, 7 = 1) ~ ^ log 2 [L(N - L)/N}. 



(8) 



Moreover, the dependence of the entropy with the mag- 
netic field is given by 

1 



S LtN (h, 7=1)- S L , N (h = 0, 7 = 1) ~ - log 2 (1 



h 2 ), 
(9) 

and thus diverges, at fixed L and N, in the limit ft-* 1 . 

Let us now discuss the more general situation 7^1, for 
which no simple analytical solution exists. In this case, 
the ground state is a superposition of Dicke states with 



3 



coefficients that can be easily determined by numerical 
diagonalizations. Upon tracing out (N — L) spins, each 
Dicke state decomposes as in Eq. @. It is then easy to 
build the (L + 1) x (L + 1) ground-state reduced density 
matrix and to compute its associated entropy. 
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FIG. 3: Entanglement entropy as a function of h near the 
critical point for 7 = 0. The solid line corresponds to the 
fitting law I I (It with a — 1/6. 



FIG. 2: Entanglement entropy at 7 = as a function of h for 
different values of N and L. Outside of the critical region, 
the entropy only depends on the ratio L/N. 

We have displayed in Fig. |3 the behavior of the en- 
tropy as a function of h, for different values of the ratio 
L/N and for 7 = 0. For ft. ^ 1, the entropy only de- 
pends on the ratio L/N. For any 7, at fixed L/N and 
in the limit h — > 00, the entropy goes to zero since the 
ground state becomes then fully polarized in the field di- 
rection. Note that the entropy also vanishes, at h > 1, 
in the limit L/N — ► where the entanglement properties 
are trivial. In the zero held limit, the entropy goes to 
a constant which depends on 7 and equals 1 at 7 — 
since the ground state is then a GHZ-like state made up 
of spins pointing in ±x directions. Close to criticality, 
the entropy displays a logarithmic divergence, which we 
numerically found to obey the law 



SL,N(h,~f) ~ -alog 2 |1 - h\, 



(10) 



where a is very close to 1 /6 for N, L 3> 1 as can be 
seen in Fig. [21 At the critical point, the entropy has a 
nontrivial behavior that we have studied focusing on the 
point 7 = which is representative of the class 7 / 1. 
There, the entropy also scales logarithmically with L as 
in the isotropic case, but with a different prefactor. More 
precisely, we find 

S L , N (h = 1,7 ¥= 1) ~ blog 2 [L(N - L)/N\. (11) 

For the finite-size systems investigated here, the prefactor 
varies when either the ratio L/N or 7 is changed, as can 
be seen in Fig. However, in the thermodynamical 
limit N, L > 1 (and finite L/N), it is likely that 6=1/3 
and does not depend on 7. In addition, at fixed L and 
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FIG. 4: Entanglement entropy as a function of L at the critical 
point for different 7 and N = 2000. The solid line corresponds 
to the fitting law 1 1 It with b = 1/3. 



N, the entropy also depends on the anisotropy parameter 
logarithmically as 

S L . N (h = 1,7) " S L , N (h = 1,7 = 0) ~ dlog 2 (l -7), 

(12) 

for all —1 < 7 < 1, as can be seen in Fig. |SJ Here again, 
it is likely that, in the thermodynamical limit, d has a 
simple (rational) value which, from our data, seems to be 
1/6. It is important to keep in mind that the limit 7 — > 1 
and the thermodynamical limit do not commute, so that 
Eq. (fEl|l is only valid for 7 ^ 1. 

Let us now compare these results with those found in 
the ID XY model. As for the LMG model, the XY chain 
has two different universality classes depending on the 
anisotropy parameter. At the critical point, the entropy 
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FIG. 5: Entanglement entropy at the critical point h = 1 as 
a function of 7. The solid line corresponds to the fitting law 
itHl with d= 1/6. 



has been found to behave as 0,0, 



L)/N], 



(13) 



where c is the central charge of the corresponding con- 
formal theory 30]. For the isotropic case, the critical 
model is indeed described by a free boson theory with 
cxx = 1 whereas the anisotropic case corresponds to a 
free fcrmion theory with with cxy = 1/2. It is striking 
to see that the entropy in the LMG model has the same 
logarithmic dependence with some prefactor which, as in 
the ID case, depends only on the universality classes [see 
Eqs. iJSJ and iJHJ]. Concerning the dependence with the 
magnetic field and with the anisotropy parameter, it is 
also worth noting that logarithmic behaviors ©, (|10|l . 
and l|12H are similar to those found in the XX and XY 
chain |3jJ except that the prefactor in the LMG model 
are different. On top, the behavior of this model with 
respect to majorization for 7^1 and as h departs from 
its critical value is completely analogous to the ID XY 
model namely, all the eigenvalues of the reduced 
density matrix obey a strict majorization relation as h 
grows, while for decreasing h one of these eigenvalues 
drives the system towards a GHZ-like state in such a way 
that majorization is strictly obeyed in the thermodynam- 
ical limit. This behavior implies a very strong "sense of 
order" in the ground state, in complete analogy to the 
XY chain. 

The logarithmic scaling with L of the ground-state en- 
tanglement entropy we have found for the LMG model, 
with well-defined values of the scaling coefficients in the 
whole parameter space, entices the search for a precise 
construction of an underlying one-dimensional local field 
theory. In order to further clarify this statement, we 
have explicitly analyzed the distribution and degeneracy 
of eigenvalues of the reduced density matrix of the sys- 



tem, that is, the structure of the partition function of 
the model. We have observed that, at 7 = 1, the spec- 
trum is doubly-degenerate with the exception of the first 
eigenvalue and follows a Gaussian behavior (as expected 
from the previous analytical calculations), while at 7 = 
the spectrum is doubly degenerate and equally spaced 
in logarithmic scale. The behavior of the spectrum for 
7 = accommodates to the typical structure imposed by 
the Virasoro algebra over the highest-weight operators in 
conformal field theories 0|. This observation might be 
a hint of a possible conformal invariance underlying the 
LMG model in some regions of its parameter space. 

Finally, let us mention that during the completion 
of this work, the entanglement entropy has also been 
calculated for the antiferromagnetic LMG model 01 for 
which the ground state is known exactly [2H 0] . 
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